Simulation of non-Abelian gauge theories with optical lattices 
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Many phenomena occurring in strongly corre- 
lated quantum systems still await conclusive ex- 
planations. The absence of isolated free quarks 
in nature is an example. It is attributed to quark 
confinement, whose origin is not yet understood. 
The phase diagram for nuclear matter at general 
temperatures and densities, studied in heavy-ion 
collisions, is largely conjectural. Finally, we have 
no definitive theory of high-temperature super- 
conductivity. Though we have theories that could 
underlie such physics, we lack the tools to de- 
termine the experimental consequences of these 
theories. Quantum simulators may provide such 
tools. Here we show how to engineer quantum 
simulators of non-Abelian lattice gauge theories. 
The systems we consider have several applica- 
tions: they can be used to mimic quark confine- 
ment or to study dimer and valence-bond states 
(which may be relevant for high-temperature su- 
perconductors) . 

Gauge theories (GT) pQ provide the basis of modern 
physics. In the "Standard Model of Particle Physics" [2] 
GT describe three of the four fundamental interactions 
(namely "electromagnetic", "weak", and "strong" inter- 
actions. The last is described by the GT known as Quan- 
tum Chromo-Dynamics (QCD) [3]). Gauge symmetry 
plays also a role in General Relativity. At the same time, 
GT are present in many effective models of condensed 
matter, e.g. antiferromagnets [3] and high-temperature 
superconductors [5J H] . Recently, the study of phase dia- 
grams of various GT has gained new attention, because of 
the discovery of topological order. Due to their stability 
against perturbations, topologically-ordered phases may 
lead to feasible quantum computers [7HTT]. 

Despite the enormous importance of GT, they defy so- 
lution. Wilson's formulation [12] of lattice gauge theo- 
ries (LGT), where continuous space-time is replaced by 
a discrete set of points, provided the first computational 
tool to study the strong coupling regime. Monte-Carlo 
(MC) simulations of LGT is the main tool to compare 
strongly-coupled aspects of QCD with experiments [15] . 
What is hard or impossible to compute with MC remains 
out of reach. For example, the mechanism of charge con- 
finement [13] , invented to explain the absence of isolated 
quarks [T5], is still debated four decades since first pro- 
posed. Furthermore, MC simulations cannot yet pro- 
vide definite predictions for hot and dense nuclear mat- 
ter [Tni [T7] , probed by heavy nuclei collisions at CERN 



and RHIC [HI [19] . GT are also invoked in explanations 
of spin-liquid phases of antiferromagnets |3] and high- 
temperature superconductivity [20] . 

Recent progress in the experimental control of quan- 
tum systems makes possible to engineer systems that 
perfectly mimic theoretical models. This is the idea of 
quantum simulators [2LH26], whose ultimate goal is to 
simulate GT, e.g. QCD, providing experimental access 
to their phase diagrams at finite temperature and den- 
sity. A more modest goal is to emulate QCD, i.e. to find 
a model sharing its interesting properties, whose realiza- 
tion may be simpler than the full theory. The first steps 
of this emulation program were to describe quantum sim- 
ulations of Abelian LGT [2"TH3"1"] . The presence of many- 
body interactions, beyond nearest neighbors, has been 
the main technical obstacle. This obstacle has been ad- 
dressed in 28J, by using mesoscopic Rydberg gates [32] . 

Here we show how to simulate, using Rydberg atoms 
[33] [34], non-Abelian gauge magnets or link models, intro- 
duced in [3"51 - |3"T] . The models have two regimes, namely 
the strong- and weak-coupling regime. We discuss here 
the origin of charge confinement in both regimes, stress- 
ing the different physical origin in each. We propose 
how to identify the flux-tubes connecting static exter- 
nal charges in each of these regimes, and provide the 
experimental protocol to observe these flux tubes. We 
conclude by discussing a qualitative technique, based on 
energy landscapes around static charges, to identify in a 
generic LGT whether chromo-electric strings, i.e. charge 
confinement, is present. 

THE MODEL 

The quantum degrees of freedom of SU(2) gauge mag- 
nets are four states assigned to each link of the lattice. 
To have a GT, we need two sets of operators, S = { S, 3 } 
and D — { T a }, for a = 0, . . . , 4 defined by 

S(S) = |0)(0| ® exp(ia • a) + |1)(1| ® I, (1) 
5(g) = |0)(0| ® I + |1)(1| ® exp(za • a), (2) 

where a — (oji, «2> CK3) is real three-component vector 
and a — (ci, 02, 03) are the Pauli matrices, and 

F° = en® J, F 7 = o 2 ®g\ j = 1,2,3, T 5 = a 3 <E)I. (3) 

The operators in the set S produce gauge transforma- 
tions. Specifically in the absence of external charges, we 
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Figure 1. The basics of non-Abelian gauge magnets a) We 
assign to each link of an oriented lattice (left), the Hilbert 
space of two qubits, each of them represented by a sphere 
(right), b) A link may be labeled in either of two equivalent 
ways: by adding a subscript referring to the adjacent site (left 
panel), or to the adjacent plaquette (right panel). We num- 
ber these labels counterclockwise. In both cases, links possess 
the original orientation of the lattice, c) On each link, one of 
the two qubits, the "position" qubit, determines the location 
of the other, the "spin" qubit, on the link. In this way, a 
red sphere on the left of a link represents the state |0) \v), 
while a red sphere on the right of a link represents the state 
1) \v). Here |0) and |1) are eigenvectors of 03, with eigen- 
values — 1 and 1, respectively. Gauge invariance forces qubits 
adjacent to a site to form singlets (represented by yellow ovals 
around these sites). For an isolated plaquette (right panel), 
the gauge-invariant Hilbert space is two dimensional. The two 
configurations of singlets are orthogonal, because they differ 
by the value of their position qubits. d) For a generic lattice, 
the gauge condition acts on the "crosses" of links adjacent to 
a given site. This condition forces the qubits adjacent to that 
site into a singlet configuration. In the left panel, we show 
one of the two linearly-independent configurations with four 
qubits adjacent to the site making a singlet. On the right, 
we show a configuration with only two adjacent qubits to the 
site forming a singlet. 

define a set of local symmetry operators 

G s (a) = 3(a) Sl ® S(a) S2 <g> S(c?) S3 <g) 3(a) S4 , (4) 

on the "crosses", which consist of the four spin states 
adjacent to given site s (labeled si , . . . , S4. See Fig. [l]b)). 
The gauge-invariance constraint selects physical states 

{ G s (a) |V) = |V) , Vs, a. (5) 

It is sufficient to impose this condition for a equal to 
x = (1, 0, 0), y = (0, 1, 0), and z = (0, 0, 1). 

The Hamiltonian H must commute with all the sym- 
metry constraints Q. As explained in detail in |36j . H is 



obtained by introducing parallel-transport operators on 
each link. The space of operators D is extended to in- 
clude operators U, which act the auxiliary spin- 1/2 space 
V. Explicitly 

3 

U = T a <E> 1 - ® r j , (6) 

3=1 

Here t j , j = 1, 2, 3 is a Pauli matrix on V. The commu- 
tation relations of the operators U and 3, 3 are 

[3(a), U] = exp(ia-T)U, [U, 3(a)] = £/exp(— ia-f), (7) 

as expected for a consistent non-Abelian GT. We now 
have enough machinery to write the Hamiltonian: 

f = J2 r i+-J2 tT ( u pi® u pi® u P3 f ® u p^+ H - c -> ( 8 ) 
1 9 P 

where A is an energy scale, I denotes the links of the 
lattice, and p denotes the elementary plaquettes of the 
lattice (pi , . . . , P4 being the links around the plaquette 
p, ordered as in Fig. [l]b)). The trace refers to the V- 
space. The coupling constant g determines if the system 
is in the weak- or in the strong- coupling regime {g — > 
and g — ¥ 00, respectively). 

In the presence of a static charge at the site s, the 
symmetry constraint at that site is modified to 

G s (d) |V) = exp (id ■ S) |V) , Va, (9) 

where S is the spin-representation of possible charges, 
1/2,1,3/2,2. Here we focus on the case of spin-1/2 
charges, i.e. S = a. In this case the charge is encoded 
by a qubit located at s and ^ can be expressed as 

Gr*(a)|V) = |V),Va (10) 

with Gf *(a) = Gs(d) ® exp (id ■ B s ). 

THE CONFINEMENT PHASE 

The form of operators S and D suggests writing the 
Hilbert space of one link &sV<S>S, which is isomorphic to 
C 2 <E)C 2 . We interpret this decomposition as describing 
qubits (the right factor, S) moving between the two ends 
of the links (the left factor, V). In accord with equations 
0, ((2), and ([§, we identify the basis of V, |0) and |1), 
with the left-end (lower-end) or the right-end (upper-end) 
of a link in the x (y) direction. In this way, we represent 
the qubit |0) \v) (|1) \v)) as a solid dot on the left (right) 
part of the link, cf. Fig[l]c). The operator 3(a) acts 
on those vectors in the left (down) two-dimensional sub- 
space |0) |i>) of the x (y) oriented link (rotating them by 
exp (id -a)). The operator 3(a) acts similarly on the 
other subspace of that link. Hence, the physical-state 
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condition ([5| forces the total spin of the qubits adjacent 
to the site s to be zero, i.e. the creation of singlet among 
pairs of those, = ^75 (tl - It), cf. Fig. [l]d). 

At weak-coupling, H in ^ reduces to its plaquette 
terms. In analogy with the Abelian case [2SJ, we exploit 
the bi-partite nature of the lattice. We imagine coloring 
the plaquettes red and black in a checker-board pattern. 
Next we consider the Hamiltonian (|Sj) , but including only 
half the terms, e.g. those on the black plaquettes. With 
this choice, the model becomes exactly solvable. As il- 
lustrated in Fig. [2]a), the ground state |O} is a product 
state of single plaquette configurations, 

i^>o = n \<f>) P > with = (iv + \pp)) . ( n ) 



id 



l°) P2 10) ps 



l)j>4 ^Pl,P2^P3,P4> anC ^ 



I 1 ),,, \°) V s S Pi,Pi S P2,P 3 - We separate the 



where \p p ) 

\K) = i-/ P2 i-/p 3 l"/p 4 

"position" part of the Hilbert space from its "spin" part 

by writing the states as elements of T 3 ® 4 ® S® A . Both \p) 

and I A) are represented in Fig. [l]c). The state |fl) fac- 

torizes into resonating-dimer states. Furthermore, each 

link, as a consequence of gauge invariance, is entangled 

with the rest of the system. 

We now turn to confinement. Adding a pair of spa- 
tially separated external charges rearranges the singlets 
into strings connecting the charges. Each string causes 
long-range entanglement (LRE) between the charges, a 
distinctive feature of non-Abelian LGT. In an Abelian 
LGT, indeed, a single string does not induce entangle- 
ment, since typically it involves flipping a line of spins 
[35]. There, the unique source of LRE between charges 
is caused by the linear superposition of several orthogonal 
string states, present also here. 

The ground state is indeed a superposition of string 
states along paths determined by both gauge invariance 
and energy minimization. A string passing through a pla- 
quette increases its energy by 5E cx A/g. Hence, strings 
touch as few plaquettes as possible. The number of ex- 
cited plaquettes is proportional to the inter-charge dis- 
tance, i.e. the charges are confined with a string tension 
proportional to A/g . This phenomenon is equivalent to 
the chromo-electric flux-tube expected in QCD between 
two colored charges. The simplest system exhibiting such 
behavior consists of only two plaquettes (Fig. |2|b). 

At strong- coupling, the plaquette term in pF may be 
neglected. The ground state is the configuration with all 
the position qubits in the state |0). Hence at any site s 
the spin qubits on S3, S4 form a singlet, see Fig[2]c) left. 
This is a product state of entangled "half-plaquettes" , of 
the form n s |0) S3 |0) S4 5 53lS4 . 

If we now insert two static charges, a line of singlets 
must readjust. As consequence of (10), the two spin 
qubits Qs 3 , Qs 4 i originally forming the singlet at s, rear- 
range. One qubit, say Qs i forms a singlet with the ex- 
ternal charge, while the position state associated to the 




Figure 2. Confinement of charges, a), b) Weak coupling, 
a) In absence of external charges, dimers resonate on configu- 
rations allowed by gauge invariance. We show two such dimer 
configurations on the left side. The lowest-energy eigenstate 
is a product state of plaquette made of resonating dimers, 
depicted as yellow spheres, b) When two static charges are 
inserted, plaquettes are entangled by strings of singlets. These 
strings resonate among gauge-invariant minimal-energy con- 
figurations (cf . left panel) . In this way, a "macromolecule" or 
"polymer" , as large as the separation of the charges, is formed. 
We show the polymer in the right panel as a yellow oval. The 
energy of such a state rises linearly with the inter-charge sepa- 
ration, thereby confining charges, c) Strong-coupling. The 
ground state of Q is the ordered configuration where all posi- 
tion qubits are in the state |0) and the spin qubits are forced 
by |5| to form singlets. This can be visualized by covering 
the lattice with singlets each formed by the two qubits posi- 
tioned left-down on the x-y links exiting a site (left panel). 
When two external charges are inserted (black spheres), the 
singlets have to rearrange in order to respect the gauge con- 
straint (10 1 (right panel). Some of the position qubits are 
flipped to |1) with an associated energy cost. Pictorially the 
corresponding spin qubits shift right-up on the x-y links. A 
string can be identified by position qubits in state jl), that is 
those spin qubits placed on the right-up of x-y links (a shaded 
blue line). The energy cost of a string is proportional to its 
length, hence, the charges are confined. 



other qubit Qs 3 changes from |0) to |1), i.e. the qubit 
moves to the opposite end of the link. There, by jBJ, 
it is forced to form a singlet with one of the two qubits 
of the same cross. The process repeats until one of the 
displaced qubits reaches the the second external charge. 
The result is a gauge-invariant string, stretching between 
the charges. The string is a line of qubits in state |1). 
The energy of each of them increases by 2A; therefore, 
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the static charges are confined. Note that the the ac- 
tual ground state is a superposition of orthogonal strings 
states inducing LRE between the external charges. How- 
ever, each single string already entangles the two charges, 
footprint of a non-Abelian LGT. 



REALIZATION THROUGH RYDBERG ATOMS 

To properly define the physical-state condition Q and 
the plaquette terms in the Hamiltonian d8|), we require 
eight-qubit interactions on neighboring sites. We can ob- 
tain long-range dipole interactions in an optical lattice by 
exciting atoms to Rydberg states. With single-site ad- 
dressing, the Rydberg dipole can be switched on and off 
at will. This property is basic to the mesoscopic Rydberg 
gate of [331 [33] . This gate can conditionally transfer the 
population between two atomic levels of a multi-atomic 
system at once, depending on the state of a "control" 
atom. Here we exploit this gate to engineer both a digi- 
tal and a mixed digital-analog quantum simulator of the 
Hamiltonian p|. 

The qubits are atoms with two (meta)stable ground 
states trapped on the links of an holographically designed 
optical lattice [32], Fig [3] a). Each link contains two 
atoms (red spheres representing both the position and 
the spin qubits). Additional two- level atoms are placed 
at the sites and at the center of the plaquettes of the 
lattice (blue spheres in Fig. [3] a)) to control the Rydberg 
gates. Geometrically, eight ensemble atoms, belonging 
to the four links around a plaquette or entering a site, 
are positioned inside the Rydberg blockade of the con- 
trols (shaded cyan and orange regions in Fig. [3] a)). All 
the atoms are in a Mott state, and are thereby fixed in 
position. 

The experiment we propose (similar to the one in [28] ) 
requires three phases: one can first prepare the ground 
state of the system for any g, both in i) absence and 
in ii) presence of static external charges and finally iii) 
observe the confinement of external charges by unveiling 
the chromo-electric strings joining them and measuring 
their tension. 

At weak-coupling, for the model with half of the pla- 
quettes on, i) is achieved by adiabatic evolution. The 
starting configuration can be chosen, e.g. as the prod- 
uct state IIi h (l+) ® l 1 ))IL„(l+) ® where l h and 
l v are the horizontal and the vertical links of the lat- 
tice, respectively, and |+) = ^(|0) + |1)), so that it 
has non-zero overlap with |fio) of (11). The evolution 



is performed for a time T (oc g/ A) using the Hamil- 
tonian H(t) = (1 - \{t))H a + X(t)(H + H G ), with the 
smooth function A(i) that fulfills A(0) = 0,A(T) = 1, 
and H a is a gapped Hamiltonian having as a unique 
ground state the starting configuration. The presence 
of H G = -AJ2 s ,a=x,yj( G s(®) + G s (d) r ) enforces gauge 





Figure 3. Lattices of Rydberg atoms needed for simulating 
SU(2)-gauge magnets, a) Holographic lattice scheme for the 
simulation of the pure LGT at weak coupling. Blue spheres 
represent the "control" two-level atoms, while red spheres are 
the "ensemble" two-level atoms. Each control atom is sur- 
rounded by eight ensemble atoms. The latter are arranged 
geometrically inside the blockade radius of the control atom. 
The blockade radii for a plaquette-control atom and the cross- 
control atom are shaded in cyan and orange, respectively. 
Single-site addressing is essential for acting selectively on each 
atom, b) An different lattice scheme has to be used to intro- 
duce static charges. It is made from super-sites including 
four two-level atoms (shaded green circle). One of these is 
a cross-control atom (blue) while the other three are used to 
encode static charges (black). Single-site addressing makes it 
possible to use i) only one of the charge qubits (having mat- 
ter with spin 1/2) or ii) several charge qubits (representing 
static charges with possible total spin 1/2, 1,3/2). The links 
of the gauge magnets, which connect these effective sites, have 
two atoms each (red spheres on shaded green oval). Inside 
the blockade radius of the cross-control atom, there are now 
eight atoms encoding the gauge magnet constituents (on the 
four links entering the site) and three charge qubits at the 
super-site (the region shaded in orange). Inside the block- 
ade radius of the plaquette-control atom, there are several 
Rydberg atoms. Using single-site addressing, some of them 
are not used by the mesoscopic Rydberg gates needed to en- 
code plaquette terms of H in Q. Indeed we selectively act 
only on the eight atoms of the gauge-magnet links around 
that plaquette (the blockade region is shaded in cyan), c) 
At strong coupling, the lattice can be considerably simplified. 
The plaquette-control atoms are absent (the Hamiltonian has 
no plaquette term) so that the number of control atoms is 
halved. We can also half the number of atoms per link, re- 
ducing it to one atom per link. Indeed, we replace the atoms 
encoding the position qubits with a double-well potential at 
each link; following (JsJ> the relative height of the two potential 
wells differs by 2A. Since spin qubits in the ground state are 
in |0), corresponding to the left (lower) well of the potential 
for x (y) links, the right (up) well is empty and the system 
is at half-filling. In order to implement (|5| the blockade ra- 
dius of the cross-controls can be limited to the first four wells 
around a site (shaded in cyan). 
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invariance, A ~ A/g. 

The adiabatic approximation is justified by choosing 
T large enough compared to the inverse of the gaps of 
H and Hq, which are both finite. The time-evolution 
operator acting on a state is implemented by a sequence 
of basic eight-qubit controlled rotations. Each of them 
involves two Rydberg gates and several single-qubit ro- 
tations (cf. [25]). The number of terms in this sequence 
poses an upper bound on the temperature of the system 
(see Methods). 

The same procedure is used to obtain ii), the ground 
state with two, or more, static charges. Their presence is 
accounted by placing new atoms on the sites of the lattice 



and enforcing ( 10 1. A possible lattice scheme that accom- 



modates several choices of static charges is presented in 
Fig. [3]b). There the sites are made by four neighboring 
potential wells and can host, apart from the control, up 
to three further atoms. The lattice spacing is chosen such 
that the desired atoms are inside the blockade radius of 
the respective controls (shaded in Fig. |3]b) ). 

Regarding iii) , we identify two strategies to detect con- 
finement. At a qualitative level, the string of singlets 
joining the two charges can be detected by using spin- 
polarization spectroscopy |40) . At a quantitative level, 
a direct measurement of the energy can be performed 
following the same approach used to perform the time 
evolution and reading out by fluorescence [H] the con- 
trol qubits. 

What we have described so far is experimentally chal- 
lenging but can be used to probe confinement at g — > 
with all the plaquettes on. This has been conjectured 
in [36] . Experimentally confinement would manifest 
through the shape of the energy-excess around the static 
charges, as an energy concentration along a tube join- 
ing the two charges. On the contrary, free charges would 
lead to spreading of the energy, footprint of a dipolar 
potential. 

In the strongly- coupled regime, we can perform a sim- 
plified experiment. Energy depends only on the position 
qubit of each link, with the state |0) (|1)) favored (penal- 
ized) by A. Hence, we perform a partial analog encoding 
by placing the spin qubits in a super-lattice producing a 
double-well potential V on each link, see Fig. |3]c). The 
two wells play the role of the position qubits and have an 
off-set of 2A. 

When two static charges are added to the system, we 
need to introduce two extra atoms above the half filled 
ground state. The creation of the strings described in 
the previous section can be mimicked by driving the sys- 
tem with e~ l $° x (t)H G dt^ £± > 2A, while simultaneously 
inducing in-well atom oscillations via AC-shaking [32J of 



the lattice at 45°. Adjusting the intensity of the shaking, 
we allow for the adiabatic adjustment of the atoms, which 
then freeze in minimal energy configurations compatible 
with gauge invariance. The strings can be observed by di- 
rect imaging of the atoms positions, e.g., by joining those 
atoms found at the right (up) end of the x (y) links. 

Note that the number of sequential Rydberg gates to 
be performed in a time step of the state preparation is 
36, i.e., small enough to satisfy A 3> thermal excitations 
(~ 1/100 Rydberg gate times) in current experiments 
(assuming a Rydberg gate frequency ~ 1/^s, and an ex- 
perimental temperature ~ 10~ 8 °K). 

Summarizing, we have proposed a quantum simulator 
of non-Abelian LGT, based on non-Abelian gauge mag- 
nets (for similar proposals see [331111]). We have identi- 
fied the mechanism producing charge confinement in this 
model. We have characterized its intrinsic non-Abelian 
nature through the part of the LRE generated by a single 
string configuration. We have designed experiments that 
allow to prepare the ground state and mimic the physics 
of charge confinement. This is only the first step towards 
the full quantum simulation of full-fledged QCD. Note 
that a slight modification of the model considered here 
allows for a relativistic dispersion relation as required by 
QCD [36]. The physics we described here is dominated 
by the presence of singlets, that play a fundamental role 
in high-temperature superconductivity. We foresee that 
the experiments we propose provide also new insights in 
this area (see also [45l [46] ) . An interesting development 
would be to apply the ideas of [47] [48] to measure in 
the experiments the intrinsic LRE carried by a single 
chromo-electric string. 

We acknowledge support from TOQATA (FIS2008- 
00784), FP7-PEOPLE-2010-IIF ENGAGES 273524, 
ERC QUAGATUA, and EU AQUTE. 



METHODS 
Quantum gates for Coherent evolution 

In the following, we give the explicit form of the driv- 
ing W and its engineering in terms of Rydberg gates. 
In particular, we discuss the upper bound on the energy 
scale of the system required by the validity of the Trotter 
approximation. As the consequence, such bound trans- 
lates in upper bound on the temperature of the atomic 
sample. 

Let us start with plaquctte term. From the explicit 
form of the operator U = oi ® a ® r° — a 2 ® 

cr-? ® , it follows that the Hamiltonian, for g —> 0, is 
H = — Y] H v , with 



H p = tr(U <g> U ® £/ f ® C/ f 



G 



cr ® cr J + 5perm. 



- ^ ( CT 1 ® cr 2 <g ' 3 <X> CT° ® CT* (g) cr J ' ' 



a* 9 ® a* - £(cr* ® c^')® 2 

) cr fc + CT 2 ® CTl ® CT2 1212 ® cr* ® cr° ® cr J ® cr fc ) + 



) CT° ® a k + <7 2 ® 3 ® 0"! ® CT* ® cr J (g) cr fc (g) cr°) 



(12) 



where cr° is the 2x2 identity matrix and e ljk is the Levi- 
Civita symbol. The indices i, j, k in the sums go from 1 
to 3. For simplicity, in the above expression we omit sub- 
scripts indicating links and plaquette. For convenience, 
the operator H p is written in the form "P® 4 ® iS® 4 . 



Hence, H p is a sum of 64 monomial tensor products 
of er^ 1 , /i = 0,1,2,3, that are not commuting, although 



some pairs of them are. Formally, we can write 



64 
,7=1 



(13) 



Let us now consider the gauge condition |5|. Defining 

H G = -AE 8 ff s = -AE., 3 =ai,«« 1 ai(G f .(3)+G.(3) t ) ) 
we have an adjoint operator by construction that can 
serve as part of the Hamiltonian. It follows that also H s 
can be written as a sum of monomials of a^. Indeed, we 
have 



1 + cos a 1 — cos a 
^ CT o H 



03 ® cr + isina- 



&0 — &3 



1 + cos a 1 — cos a 
s CT o ^ 



03 ® c + i sin a 



crM +ff.C. (14) 



As the decomposition is a very lengthy expression, we 
count the total number of monomials. The Hermi- 
tian terms are the ones containing even powers of sin a. 
For a given j, we have (q)2 4 x (sin a) terms, plus 
(!)2 4 x (sin a) 2 terms, plus (|)2 4 x (sin a) 4 terms. As 
the (sin a)°-terms are the same for any j = 1,2,3, the 
total number of terms is (1 + 3 x 7)2 4 = 352. Hence, 
formally we can write 

352 8 

^£q1 j) , = (15) 

J=l n=l 

As explained in the main text, the adiabatic prepa- 
ration requires also an auxiliary gapped Hamiltonian H a 
that has the configuration we start with as unique ground 
state. Such a H a ensures that the unitary evolution 

t=T 

t=o 



t=o 

is adiabatic for X(t) -C 1, i.e. i « T. A suitable choice 
for the initial state is I 1 )) Hi v (|+) ® |0)), where 

lh and l v are the horizontal and the vertical links of the 
lattice, respectively, and |+) = r^(|0) + |1)), so that it 
has non-zero overlap with |f2 ) °f ( |H[ ) ■ A consequent 
suitable choice for H a is 

H a = —A I £((7i ® CT° - cr «) cr 1 ) £(oi ® cr° - CT ® cr 1 ) 

which is very easy to implement in terms of single-qubit 
rotations. 

It is worth to notice that alternative plaquette terms 
(to the one we consider here) and the electric term 
A ^ s T 5 = J^i a 3 ® a ° can be also easily decomposed in 
monomials of cr M and implemented as explained below. 
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Now we are able to compute the driving operator W 
in terms of the Q's. Indeed, we have 



if t t+H (i-X(t'))H a dt' 



W{t) 



-ij t t+it \(t')AJ2 s 



dt' -i f. 



as the dynamics and, partially, the structure of Hilbert 
space are encoded directly in the lattice, we have to per- 
/ form only the coherent evolution of an exponentially sim- 
pler gauge condition: the cross operator acts on only four 
' E p ,j | 3p J)(lt ' i qubits as in the Abelian case [55], instead of eight. 
(lTfhe new H s is 



-if*+" A(i')(A£ a H„ + f E p H p )dt' 



since the two Hamiltonians commutes by construction, 
due to ([7]). The driving is efficient as, for A, A =^ 0, H 
and Hq are both gapped, with gaps proportional to A 
and A, respectively. 

Note that the scale A has to chosen appropriately. In 
absence of charges, as \Hq, H] = and |f2o) is the simul- 
taneous minimum of both H and Hq, it is sufficient to 
take A ~ A/g. However, when charges are present, the 
ground state is an eigenvector of iJ, but not with min- 
imal eigenvalue, if we do not limit the Hilbert space to 
the gauge-invariant subspace. We ensure that this is the 
case by requiring that the ground state minimizes Hq, 
i.e. requiring A sufficiently greater than A/5, precisely 
A>||#||. 

As explained in [33J and [55J, any term of the form 
e iQ< 2 ( ' can be implemented, up to single-qubit rotations, 
by applying the mesoscopic Rydberg gate twice. Hence, 
we can expand in Trotter approximation any of the two 



exponentials in (17 1 in product of exp of the Q's and 



realize the evolution as a sequence of gates. For instance, 



no 



A y 



N 



N 



for -~— < 1. 
N q 



(18) 



The crucial point to be noticed is that N cannot be ar- 
bitrary large. In order to perform a time-step of evo- 
lution, our Rydberg apparatus physically needs a time 
T step = 2# of Q's xtfj, where ta is the delay of a single 
Rydberg gate. It follows that 



1 



T, 



step 



(19) 



which for the plaquette term above, without further opti- 
■ ~ 10 2 KHz, while for the gauge 



mization, is of order 
term is about five time smaller. This means that ground 
state and dynamics in a generic regime of (|| requires 
a considerably cold atomic sample, about one order of 
magnitude colder than in present experiments. However, 
we do not discard that the above implementation can be 
improved and made less demanding applying gate opti- 
mization. 



Strong coupling experiment 

The situation improves considerably in the simplified 
scheme we propose in the strong coupling limit. Indeed, 



H s = J2( eia3 ) m = 2(cosa) 4 a 0$ 



— 2(cosasina) 2 (cr : '^'" ® ^o® 2 _j_ 5p erm ) _|_ ( s i n o!) 4 <7 

and the only non-trivial terms that we have to engi- 
neer with Rydberg gates are 18, i.e. T t 



4 ,84 



step 



36t ff . As 



briefly mentioned in the main text, such T step is suffi- 
ciently smaller than the typical time-scale associated to 
(the inverse of) thermal fluctuations in current state-of- 
art experiments, T < 10~ 8 °K. 

We conclude with few comments on AC-shaking pro- 
cedure. As explained in the main text, we allow the sys- 
tem to reach the minimal energy configuration during the 
adiabatic evolution by inducing an hopping between the 
wells forming the links. In this way, the value of the po- 
sition qubit, i.e. the position of the two-levels atom in 
each well, can change to the less energetic one while sat- 
isfying the new gauge condition in presence of charges. 
The gauge condition determines the configuration of the 
gauge qubits, forming singlets between pairs of them. An 
AC-periodic forcing of the lattice can provide this hop- 
ping in a selective way. Indeed, the AC-induced hopping 
is obtained only if the frequency of the forcing u> is com- 
mensurable with 2 A, the relative energy off-set between 
the two wells, i.e. A = %hui, v integer [33]. This is suffi- 
cient to avoid the generation of unwanted tunneling terms 
as, for instance, of the control atom to the link wells and 
vice-versa. Furthermore, other unwanted hopping terms 
between different links are strongly suppressed by the 
exponential decay of Wannier functions. In order to in- 
duce the same rate of in-well oscillations for horizontal 
and vertical links, the periodic forcing F is chosen to be 
parallel to the plane at a 45° angle, F \\ (1, 1,0). 
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